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Abstract. We test the accuracy of the ALI method, widely used in stel- 
lar atmosphere modelling, by solving exactly a standard radiative transfer 
problem in plane-parallel geometry. Some recommendations are given for a 
practical use of this method. 

1 Introduction 



We are interested in solving the radiative transfer equation in a homogeneous 
and isothermal plane-parallel atmosphere in local thermodynamical equilibrium, 
*^ assuming isotropic and monochromatic scattering. If no radiation is incident on 
the boundary layers, the thermal source function is eB(T), where e is the photon 
destruction probability per scattering and B(T) is the Planck function at temper- 
ature T. The source function S, normalized to the Planck function, satisfies the 
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^ ' following integral equation 

S(r)=e+(l-e)(AS)(r), (1.1) 
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where 



(A5)(r) = ijf E^T-T'DS^dT'. (1.2) 

E\{t) — J exp(— r/pi) dfi/fi is the first exponential integral function, and the 
optical depth variable r covers the range [0,r*], r* being the optical thickness of 
the atmosphere at a given frequency. Equation (|1.1|) models the multiple scattering 
of photons in a continuum or in a spectral line with a rectangular profile. 
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Chevallier and Rutily (2003) have solved recently Eq. (|l.lfl with an accuracy 
better than 10~ 10 for any value of e, r* and r. Their solution can safely be used 
as a reference solution for testing the accuracy of a numerical code, e.g., the ALI 
code widely used in stellar atmosphere modelling. Our focus is to present some 
tests for this code, which we first describe in the next section. Some additional 
tests can be found in Chevallier et al. (2003). 



2 Brief description of our ALI code 



The ALI code used in this paper is a combination of an accelerated iterative 
method (with a diagonal A-operator) and a formal solver based on parabolic short 
characteristics. Recent reviews on this approach are Paletou (2001) and Hubeny 
(2003). The grid we chose for the integration with respect to r over [0, r*] is 
logarithmic and symmetric with respect to the mid-plane, with n T = 9 points per 
decade starting from r = 10~ 4 and including the r = point. The A-operator is 
defined by an integration over a directional variable [i in [— 1 , 1] , which is performed 
with the help of a Gaussian quadrature (n M = 5 Gauss-Legendre points on [0, 1], 
plus opposite points). The iterations can be stopped after N > N c iterations, 
where N c is the number of iterations required for the accuracy c?m (i-e., maximal 
relative error) to reach the convergence threshold (Fig. 1). As N c is not known a 
priori, the iterations are usually stopped when the relative change R c of the source 
function between two successive iterations is less than a given small number, for 
instance 10~ 3 . 




Fig. 1. Accuracy <1m and relative change R c of the ALI code as a function of the number 
of iterations N for a "strong" spectral line e = 10 -8 , t" — 2 x 10 8 . 



3 The accuracy tests 

Figure 2 shows the relative error of our ALI code as a function of r for N = 
1000 > N c . The three curves correspond to typical values of (e, t*) in a continuum 
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(e = 0.01, t* = 20), in an "average" spectral line (e = 10~ 4 ,t* = 2000) and in 
a "strong" spectral line (e = 10 _8 ,t* = 2 x 10 s ). The relative error is maximal 
when t is smaller than the thermalization depth l/fc(e), equal to 6, 58 and 5774 
for the continuum, average line and strong line respectively (Chevallier and Rutily 
2003). It significantly improves in the deep layers of the star, up to 10 -8 . We 
confirm that the accuracy of our ALI code for N > N c does not depend on e and 
t* , but only on the refinement of the r-grid and the ^-grid. 

Figure 3 shows the number of iterations N c needed to achieve convergence, 
for an extended range of (e, r*). We note an important slowing down of the 
convergence when e — > and t* — > oo. From preceding studies of the semi-infinite 
case, it is well known that N c strongly depends on e (Trujillo Bueno & Fabiani 
Bendicho 1995). The present study of the finite case clarifies the influence of r*: 
N c depends only on e when fc(e)r* > 100 and drops down for weaker values of r*. 




Fig. 2. Relative error d(r) of the ALI code for N = 1000 > N c . Black dots show the 
thermalization depth l/k{e) for each case (k(e) ~ \/3e when e — *• 0). 

In practice, an accuracy better than 1% is generally required while choosing 
R c = 10 -3 as the iteration stopping criterion. Figure 4 shows the accuracy of ALI 
for a set of slabs such as e = 10~ p and t — 2 x 10 p (p = 1, . . . , 8) when R c = 10~ 3 
and R c — 10~ 4 . The criterion R c = 10 -3 does not warrant an accuracy of 1% for 
any p. We thus recommend the value R c < 10~ 4 to ensure an accuracy of 1% in 
stellar atmosphere calculations. 

4 Conclusion 

In this particular problem (|1.1|) . our code was able to converge in all situations, its 
accuracy being bounded by 5 x 10 -3 and 10 -2 when using the stopping criterion 
R c = 10~ 4 . We intend to extend this work to more realistic problems (e.g., non 
isothermal atmospheres), replacing the Jacobi scheme by Gauss-Seidel or Succes- 
sive Over- Relaxation schemes (see Trujillo Bueno & Fabiani Bendicho 1995). 
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Fig. 4. Accuracy dwi as a function of parameter p, where e = 10 p and r* = 2 x 10 p , for 
R c = 10~ 3 and R c = 10 -4 . The dotted- line figures out the 1% accuracy level. 
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